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ABSTRACT 

By combining Herschel-^V^RE observations obtained as part of the Herschel Virgo Cluster Survey with 21 cm Hi data from the literature, we 
investigate the role of the cluster environment on the dust content of Virgo spiral galaxies. We show for the first time that the extent of the dust disk 
is significantly reduced in Hi-deficient galaxies, following remarkably well the observed 'truncation' of the Hi disk. The ratio of the submillimetre- 
to-optical diameter correlates with the Hi-deficiency, suggesting that the cluster environment is able to strip dust as well as gas. These results 
provide important insights not only into the evolution of cluster galaxies but also into the metal enrichment of the intra-cluster medium. 
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1. Introduction 

It is now well established that the evolution of spiral galaxies 
significantly depends on the environme nt they inhabit. Th e re- 
duction in the star formation rate (e.g.. iLewis et al.ll2002h and 
atomic hydrogen (Hi) content (e.g. jGiovanelli & Haynes'l 19851) 
of galaxies when moving from low- to high-density environ- 
ments indicates that clusters are extremely hostile places for 
star- forming galaxies. However, a detailed knowledge of the ef- 
fects of the environment on all the components of the interstellar 
medium (ISM) is still lacking. Particularly important is our un- 
derstanding of how the environment is able to aflTect the dust con- 
tent of cluster spirals. Dust plays an important role in the process 
of star formation, since it acts as a catalyzer for the formation of 
molecular hydrogen (H2, from which stars are formed) and pre- 
vents its dissociation by the interstellar radiation field. Thus, the 
stripping of dust might significantly aflTect the properties of the 
ISM in inf ailing cluster spirals. 

Since dust is generally associated with the cold gas compo- 
nent of the ISM, it is expected that when the Hi is stripped part 
of the dust will be removed as well, but no definitive evidence 
of a reduced dust content in cluster galaxies has been found so 
far. For a fixed morphological type, Hi-deficient galaxie£l ap- 
pear to have higher IRAS /(1 00 /u rn)/ f (60 /urn) flu x density 
ratios (i.e., colder dust temperatures. iBicay & Giova nelli 198'^ 
and lower far-infrared (FIR) flux densities per unit optical area 
(■Doyon & Joseph. 1989) than gas-rich galaxie s. However by u s- 
ing ISO observations of the Virgo cluster (iTuflTs et al.l I2002D . 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 

^ The Hi-deficiency (defm) is defined as the diff'erence, in logarith- 
mic units, between the observed Hi mass and the value expected from an 
isolated galaxy with the sam e morphological type and optical diameter 
(iHavnes & Giovanellil[T984l) . 



iPopescu et al.l (l2002h find no strong variation with cluster- 
centric distance in the dust properties of each morpholog- 
ical type. Only the most extreme Hi-deficient galaxies ap- 
pear to be lacking a c old dust component. More recently, 
iBoselli & Gavazzil (l2006l) have revealed an interesting trend of 
decreasing dust masses per unit of H-band luminosity with de- 
creasing distance from the center of Virgo. Thus, it is still an 
open issue whether or not dust is removed from infalling cluster 
spirals. 

The launch of Herschel (iPilbratt et al.l 1201 Oh has opened a 
new era in the study of environmental eflPects on dust. Thanks 
to its high spatial resolution and sensitivity to all dust com- 
ponents, Herschel will be able to determine if cluster galaxies 
have lost a significant amount of their dust content. Ideally, this 
analysis should be done on a large, statistically complete sam- 
ple, following the same criteria used to define the Hi-deficiency 
parameter ( Hay nes & Giovanellil 1198 4*): i.e., by comparing the 
dust content of galaxies of the same morphological type but 
in diflTerent environments. By observing a significant fraction 
(-64 deg2) of the Virgo cluster at 100, 160, 250, 350 and 500 
/im, the Herschel Virgo Cluster Survey (HeViCS, IPavies et"aD 
2010, hereafter Paper I; see also http://www.hevics.org) will 
soon provide the optimal sample for such an investigation. 
In the meantime, with the first HeViCS data it is possible to 
use a more indirect approach and compare the extent of the 
dust disk in gas-rich and gas-poor cluster galaxies. Since pre- 
vious studies have shown that the Hi stripping is associated 
with a ' truncation Q of the ^as (jCav atte et al.l [l994l) and star- 
forming disk ( Koopma nn & Keniievl 2004: Ca tineUa et"ani2005l: 
iBoselli & Gavazzil i2006h . if the dust follows the atomic hydro- 



^ The term 'truncation' is used here to indicate either an abrupt steep- 
ening of the surface-brightness profile or, more simply, a significant re- 
duction in the disk scale-length compared to the optical one. 
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Fig. 1. The optical (top), 250 jim (middle) and Hi (bottom) maps for galaxies in our sample with different degrees of Hi-deficiency. 



gen we should find a reduction in the extent of the dust disk with 
increasing Hi-deficiency. 

In this paper we will take advantage of the HeViCS obser- 
vations, obtained as part of the Herschel Science Demonstration 
(SD) phase, to investigate the correlation between the dust dis- 
tribution and gas content in cluster galaxies. 

2. Observations and data reduction 

A ~245'x 230' field in the center of the Virgo clust er has 
been observed by Hersche l using the SPIRE/PACS (Griffin eran 
I2OIOI: iPoglitsch et al.ll201 0) parallel scan-map mode as part of 
the SD observations for HeViCS. In this paper we will focus 
our attention on the 3 SPIRE bands only. The full widths at half 
maximum of the SPIRE beams are 18.1, 25.2, 36.9 arcsec at 
250, 350 and 500 yum, respectively. Details about the observa- 
tions and data reduction can be found in Paper I. The typical rms 
noise across the whole image are -12, 10, 12 mJy/beam at 250, 
350 and 500 yum, respectively (i.e., ~2 times higher than the con- 
fusion noise). No spatial filtering is applied during the data re- 
duction, making SPIRE maps ideal to investigate extended sub- 
millimetre (submm) emission. The uncertainty in the flux cali- 
bration is of the order of 15%. 

In order to investigate how the dust distribution varies with 
the degree of Hi-deficiency in Virgo spirals, we restricted our 
analysis to the 15 spiral galaxies in the HeViCS SD field for 
which Hi surface density profiles are available. The Hi maps 
are obtained from the recent 'VLA Imaging of Virgo in Atomic 
gas' (VIVA) survey (Chung et al. 2009, 13 galaxies: NGC4294, 
NGC4299, NGC4330, NGC4351, NGC4380, NGC4388, 
NGC4402, NGC4424, NG C4501, NGC4567, NGC4568, 
NGC456 9, NGC4579), f rom ICavatte et al.1 ([T99l NGC4438) 
and from Warmels (198o, NGC4413). Hi-deficiencies hav e been 
determined following the prescription presented in iChung et al.l 
(I2OO9I) . This method is slightly diff'erent from the original 
definition presented by iHaynes & Giovanellil (Il984l) . as it 



assumes a mean Hi mass-diar neter relation , regardless of the 
morphological type. Following ^ Chung et al.l 12009), we use the 
diff'erence between the type-dependent and type-independent 
definitions as uncertainty in the Hi-deficiency parameter. We 
note that, on average, this value is sm aller than the intrinsi c 
scatter ofdefui for field galaxies (-0.27. lFumagalli et al.ll2009h . 

Surface brightness profiles in the three SPIRE bands were 
derived using the IRAF task ellipse. The center was fixed 
to the galaxy's optical center (taken from the NASA/IPAC 
Extragalactic DatabaseEl) and the ellipticity and position angle 
to the s ame values adopted for the Hi profiles taken from the lit - 
erature (IChung et al.ll2009l: ICavatte etal.ll 19941: IWarmelsl[T988h . 
The sky background was determined within rectangular regions 
around the galaxy and subtracted from the images before per- 
forming the ellipse fitting. Each profile was then corrected to 
the 'face-on' value using the inclinations taken from the litera- 
ture. All the galaxies in our sample are clearly resolved in all the 
three SPIRE bands: e.g., on average -4-5 beam sizes at 500 yum. 
Submm isophotal radii were determined at 6.7x10"^, 3.4x10"^ 
and 1.7x10"^ Jy arcsec"^ surface brightness level in 250, 350 
and 500 jum respectively. These are the average surface bright- 
nesses observed at the optic al radius (25 mag arcsec"^ in B band, 
Ide Vaucouleurs et anil99lh in the four non Hi-deficient galaxies 
(defui <0.3) in our sample (NGC4294, 4299, 4351, 4567) and 
correspond to ~2-3(T noise level. Of course, this choice is rather 
arbitrary and it has no real physical basis. However, as discussed 
in § 3, the result does not depend on the way in which the isopho- 
tal radii have been defined. Although many of our galaxies show 
some evidence of nuclear activity, we do not find a single case 
in which the nuclear s ubmm emission do minates the emission 
from the disk (see also ^Sauvage et al.|[2010 ). Thus, the isophotal 
radius is a fair indication of the extent of the dust disk. 



http://nedwww.ipac.caltech.edu/ 
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Fig. 2. The ratio of the submm-to-optical diameters versus Hi- 
deficiency in the three SPIRE bands. Squares are for Sa-Sab, 
stars for Sb-Sbc and hexagons for Sc and later types. For com- 
parison, the triangles show the same relation for the Hi-to-optical 
diameter ratio, where the Hi isoph otal diameters are taken at a 
surface density level of 1 M© pc"^ (IChung et al.ll2Q09h . 



3. Results & discussion 

In Fig. [T] we compare the optical, 250 yum and Hi maps for a 
subsample of our galaxies with different levels of Hi-deficiency. 
In highly deficient spirals the 250 jim emission is significantly 
less extended than the optical, following remarkably well the ob- 
served 'truncation' of the Hi dislfl. This is confirmed in Fig. [2l 
where we show the ratio of the submm-to-optical isophotal di- 
ameters as a function of defni for the 15 galaxies in our sam- 
ple. For all the three SPIRE bands we find a strong correlation 

(Spearman correlation coefficient r^ 0.87, corresponding to 

a probability P(r > r^) > 99.9% that the two variables are cor- 
related) between the submm-to-optical diameter ratio and defm- 
Although qualitatively supported by Fig.[T] this correlation alone 
does not imply a change in the shape of the submm profile. A 
decrease in the central submm surface brightness of gas-poor 
galaxies could produce a similar trend without the need to in- 
voke a reduction in the disk scale-length. However, Fig.[3]and|4l 
clearly exclude such a scenario. In Fig. [3] we show that, while the 
350 yum flux per unit of 350 yum area (i.e., the average submm sur- 
face brightness) is nearly constant across the whole sample, the 
350 //m flux per unit of optical area significantly decreases with 
increasing defni. This is even more evident in Fig. |4] where the 
average surface brightness profiles in bins of normalized radius 
for gas-rich and gas-poor galaxies {defm >0.96; i.e., NGC4380, 
NGC4388, NGC4424, NGC4438, NGC4569) are shown. While 
the central surface brightness is approximately the same, the pro- 
file of Hi-deficient galaxies is steeper than in normal galaxies and 
falls below our detection limit at approximately half the optical 
radiu£l. We can thus conclude that Hi-deficient galaxies have 
submm disks significantly less extended than the optical disks, 
following closely the 'truncation' observed in Hi. 



^ See also lPohlen et al.l (1201 Oh for an analysis of the two grand design 
Virgo spirals NGC4254 {defm ~ -0.10) and NGC4321 {defm -0.35). 

^ This also confirms that the correlation shown in Fig. [2] is not qual- 
itatively affected by the definition of submm isophotal radius adopted 
here. 
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Fig. 3. The 350 yum flux per unit of 350 yum area (left) and optical 
area (right) versus Hi-deficiency. Symbols are as in Fig.[2l 
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Fig. 4. Average submm surface-brightness profiles in bins of nor- 
malized radius for normal and highly Hi-deficient galaxies. 



Interestingly, from Fig. [2] and [3] it emerges that the extent of 
the dust disk is significantly reduced compared to the optical disk 
only for high Hi-deficiencies (J^/hi~0.8-1), i. e. when the atomic 
hydrogen starts to be stripped from inside the optical radius. 

We now need to consider whether we are just observing a 
trend due to a different mix of morphologies between gas-rich 
and gas-poor galaxies. Although Hi-deficient systems are of ear- 
lier type than gas-rich spirals, our result does not change if we 
focus our attention on Sa-Sbc galaxies only (i.e., 80% of our 
sample). Since in this range the ave rage 850 yum scale-length- 
to-optica l radius (IThomas etaP 12004) and Hi-to-optical radius 
(Cayatte etal .11 19941) ratios donot vary significantly (i.e., less 
than Icr) with galaxy type, morphology alone cannot be respon- 
sible for the correlations shown in Fig. |2]and[3] Moreover, all 
the highly Hi-deficient galaxies in our sample are well known 
perturbed Virgo spirals, on which the inff uence of the cl uster en- 
vironment has already been proven (e.g., Vollmej|2009h . So, the 
difference in the dust distribution between gas-poor and gas-rich 
spirals observed here is likely due to the eff'ect of the cluster en- 
vironment and is not just related to the intrinsic properties of 
each galaxy. Future analysis of a larger and more complete sam- 
ple will allow us to further disentangle the role of environment 
from morphology on the dust distribution in nearby spirals. 

A 'truncation' in the surface brightness profile of NGC4569 
(the most Hi-deficient galaxy in our sam ple) has alrea d y bee n 
observed at Spitzer 24 and 70 yum bv iBoseUi et al.1 (l2006h . 
However, while a reduction in the 24 and 70 yum surface bright- 
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Fig. 5. The 250 yum-to-K-band (left) and 500 yum-to-K-band flux 
density ratios versus Hi-deficiency. Symbols are as in Fig.O 



ness might just be a direct consequence of the quenching of the 
star formation in gas-poor galaxies, this scenario is not valid in 
our case. For /1~ 100-200 yum, the dust emission does not come 
predominantly from grains directly heated by photons associ- 
ated with star formation activity, but from a colder component 
heated by photons of the difl^use interstellar radiation fie ld (e.g., 
IChini et al.|[T986l: iDraine eTaDl2007l: iBendo et al.ll20TQh . Since 
this colder component dominates the dust mass budget in galax- 
ies, the trends here observed are likely not due to a reduction in 
the intensity of the ultraviolet radiation field, but they imply that 
in Hi-deficient galaxies the dust surface density in the outer parts 
of the disk is significantly lower than in normal spirals. 

An alternative way to compare the properties of normal and 
gas-poor Virgo spirals is to look at their submm-to-near-infrared 
colours. Since the K-band is an ideal proxy for the stellar mass 
and the SPIRE fluxes provide an indication of the total dust 
mass, it is interesting to investigate how the f (250)/ f(K) and 
f(500)/f(K) flux density ratios vary with defm- We find that 
highly Hi-deficient galaxies have f(250)/f(K) and f(500)/f(K) 
ratios a factor -2-3 lower than normal galaxies (Fig.O. This pro- 
vides additional support to a scenario in which gas-poor galaxies 
have also lost a significant fraction of their original dust content. 

By comparing the dust mass per unit of H-band luminosity 
for a sam ple of late-type galaxies in the Coma- AbelU 367 su- 
percluster. lContursi et al.l (2001) find no significant diff'erence in 
the dust content of normal and Hi-deficient spirals, apparently in 
contrast with our results. However, such a diff^er ence is due (at 
least i n part) to the fact that the sample used by IContursi et al.l 
(l200lh does not include any galaxy with defui >0.87. It is easy 
to see in Fig. [5] that, for defui <0.87, almost no trend is ob- 
served between the f (500)/ f(K) (or f (250)/ f(K)) flux density 
ratio and defni. In fact, the Spearman correlation coefficient 

drops from r^ 0.78 to 0.28 and 0.11 (i.e., a drop in 

the probability that two variables are correlated to 80 and 40%) 
for the f (250)/ f(K) and f (500)/ f(K) ratios, respectively. This 
implies that the two variables are no longer significantly corre- 
lated, highlighting once more that substantial dust stripping is 
observed only if the ISM is removed from well within the opti- 
cal radius. 



4. Conclusions 

In this paper, we have shown that in Hi-deficient galaxies the 
dust disk is significantly less extended than in gas-rich systems. 
This result, combined with the evidence that Hi-deficient objects 
show a reduction in their submm-to-K-band flux density ratios, 
suggests that when the atomic hydrogen is stripped part of the 
dust is removed as well. However, the dust stripping appears ef- 
ficient only when very gas-poor spirals are considered, imply- 
ing that in order to be significant the stripping has to occur well 



within the optical radius. This is consistent with iThomas et al.l 
(2004) who found that the 850 yum scale-length of nearby galax- 
ies is smaller than the Hi, suggesting that outside the optical ra- 
dius the gas-to-dust ratio is higher than in the inner parts. 

Our analysis provides evidence that the cluster environment 
is able to significantly alter the dust properties of infalling spi- 
rals. We note that this has only been possible thanks to the unique 
spatial resolution and high sensitivity in detecting cold dust pro- 
vided by the Herschel-SFIRE instrument and to the wide range 
of Hi-deficiencies covered by our sample. Once c ombined with 
the direct detection of stripped dust presented by I Cortese et al.l 
(2010) and "Gomez et al.' ("2010), our results highlight dust strip- 
ping by environmental eff'ects as an important mechanism for 
injecting dust grains into the intra-cluster medium, thus con- 
tributing to its metal enrichment. This is consistent with numer- 
ical simulations which predict that ram pressure alone can al- 
ready con^ibute_^40%of_the enrichment of the ICM in clus- 
ters (iDomainko et al.l 12006) . Interestingly, the stripped grains 
should survive in the hot ICM long en ough to be observed 
(iPopescu et al.ll2000l: [Clemens et al.ll2010l) . 

Once completed, HeViCS will allow a search for additional 
evidence of dust stripping and place important constraints on 
the amount of intra-cluster dust present in Virgo. Moreover, in 
combination with the Herschel Reference Survey (iBoselli et al.l 
2010b), it will be eventually possible to accurately quantify the 
degree of dust-deficiency in Virgo spirals. 
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